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Supernova searches have been been suggested as a method 
for determining precisely the current value and time variation 
of the equation of state, w, of the dark energy component 
responsible for the accelerated expansion of the Universe. We 
show that the method is fundamentally limited by the fact 
that luminosity distance depends on w through a multiple 
integral relation that smears out information about w and its 
time variation. The effect degrades the resolution of w that 
can be obtained from current data. 



Recent observations suggest that most of the energy 
density of the Universe consists of a dark energy com- 
ponent with negative pressure that causes the expan- 
sion rate of the Universe to accelerated A key chal- 
lenge for cosmology and for fundamental physics is to 
determine the nature of the dark energy. One possibil- 
ity is that the dark energy consists of vacuum energy 
or cosmological constant, in which case the equation of 
state is w = p/p = —1, where p is the pressure and p 
is the energy density of the dark energy. An alterna- 
tive is quintessence,0 a time-evolving, spatially inhomo- 
geneous energy component with negative pressure, such 
as a scalar field slowly rolling down a potential. For 
quintessence, the equation of state is typically a function 
of redshift, w(z), whose value differs from -1. Hence, a 
precise measurement of w today and its time variation 
could distinguish between the two possibilities and pro- 
vide important clues about the dynamical properties of 
dark energy. 

Searches for type la supernovae at deep redshift have 
provided the most direct evidence that the expansion rate 
of the Universe is accelerating.Em The supernovae appear 
to be standard candles which can be used to measure 
the luminosity distance-redshift relation. By measuring 
50 supernovae out to redshift Jiear z = 1, the Super^. 
nova Cosmology Project (SCP)B and the High -z Surveycl 
Project have each found strong evidence that the Uni- 
verse is accelerating and that the equation of state of the 
dark energy component is negative.EJ 
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FIG. 1. (a) The luminosity distance HooIl(z) for nine 
choices of equation of state wq(z) for the dark energy shown 
in (c), where Ho is the current value of the Hubble parame- 
ter. All models have Q m = 0.3. (b) Illustrates that the per- 
centage deviation of cIl(z) from a a cosmological model with 
f2 m ~ 0.3, Qq = 0.7 and wq = —0.7 = const, is less than 1%. 
If one artificially restricts wq to be constant, then the range 
of models collapses to the region between the dashed lines. 

A supernova search extended to greater z can make a 
much more precise determination of the luminosity dis- 
tance as a function of red shift Jj'Ll di,{z), perhaps to bet- 
ter than 1% uncertainty out to redshift z = 2. (1% is 
probably an optimistic estimate of the limiting system- 
atic uncertainty.) Does this enable a precise determina- 
tion of the equation of state of the dark energy compo- 
nent and its time variation? As we show in this paper, 
the answer is no. The inherent limitation is theoreti- 
cal: the luminosity distance depends on w{z) through 
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a multiple-integral relation that smears out detailed in- 
formation about w(z). Consequently, the value of w{z) 
today is poorly resolved and no useful constraint can be 
obtained about its time variation. 

The problem can be immediately appreciated from 
Fig. 1, which compares <Il(z) for an assumed cosmologi- 
cal model [Q m = 0.3, Oq = 0.7 and wq = —0.7 = const., 
where £l m ,Q is the ratio of the (matter, quintessence) en- 
ergy density to the critical density] with eight other mod- 
els chosen as examples where <Il{z) is nearly degenerate 
with the assumed model. In this figure and throughout 
the paper, we assume the Universe is cosmologically flat 
and the speed of light c = 1. (Henceforth, we use the 
subscript Q to label the dark energy component, be it 
quintessence or cosmological constant.) Figure la shows 
that cIl(z) is nearly identical for the set of models as 
individual curves can hardly be distinguished. Fig. lb 
displays the percentage deviation of <1l{z) from the as- 
sumed model, where it can be seen that the deviation 
is less than one per cent out to redshift z = 2. Fig. lc 
then shows wq(z) for the respective models. The striking 
result is the wide range of wq(z) that produces nearly 
the same di,(z) as the assumed model. If one expands 
wq(z) = Wo + w\z + W2Z 2 + . . . , then, for this particular 
collection of models, wq varies between -0.55 and -0.9 (a 
total span of 50% about the assumed value, wq = —0.7) 
and wi = diVQ/dzQ varies between -1.1 and +1.6. (The 
subscript "0" refers to present-day values of parameters.) 

Note that the degenerate models chosen for the illus- 
tration span a larger range of \dwQ/dz\ = 0(1) than 
most realistic models predict. Typically, \dwQ / dzo\ -C 1 
because wq{z) is bounded in most cases to lie between -1 
and +1 in order that the dark energy obey the positive 
energy condition and be stable under perturbations. The 
large uncertainty in w\ = dvjQ /dzQ means that little use- 
ful information is obtained about the magnitude or sign 
of the time variation of wq. Also, wq is poorly resolved. 
The resolution of wq(z) degrades significantly further if 
one includes the uncertainty in Q m , as shown in in Fig. 2 
(see discussion below) . 

Our conclusion may seem at odds with some pro- 
jections DLwhat can be obtained in future supernova 
searchesJj'Ll Many analyses assume wq — const. If we 
impose this condition, then the range of models that fit 
collapses to the narrow region between the dashed lines in 
Fig. lc, giving a misleading impression that wq(z) is well- 
resolved. However, if consideration is extended to mod- 
els in which wq is z-dependent, such as the linear form 
wq(z) = wq + W\Z, the result is dramatically different. 
A very wide range of (wo, w\) produces nearly identical 
dzy(z) because the differences are smoothed out by the 
multi integral relation derived below between luminosity 
distance and wq(z). This degeneracy accounts for the 
results found in Figs, lc and 2, but is missed if one artifi- 
cially restricts wq to be constant. Notice in Fig. lc that 
including non linear forms for wq(z) enhances the uncer- 



tainty in wq and wi even further. Among studies which 
have considered time varying wq(z), our results agree 
with sQjjuetrlilJ but seem significantly less optimistic than 
others. Qllil In the latter cases, the assumptions about the 
observations are similar but various subtle factors, such 
as the use of fitting functions rather than exact expres- 
sions for di,(z) or imposing the constraint wq > — 1, 
combine numerically to reduce artificially the degener- 
acy Extending searches to yet deeper redshift (z > 2) 
does not help either because the effect of quintessence on 
c?l(z) is proportional to Oq which becomes very small at 
deep redshift. 

The key to understanding these conclusions is the 
relation between luminosity distance and the equation 
of state. The luminosity distance is related to the 
Robertson- Walker scale factor a(t) through the equation 
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where the redshift z satisfies l + z = ao/a, and H is 
the Hubble parameter, H 2 = H$(p T (z)/p T (0)). We 
assume a flat Universe. The subscript "T" refers to 
the total equation of state of the the combined matter- 
quintessence fluid. Integrating the energy conservation 
equation, 



Pt = — 3if (1 + wt)pt 
we can reexpress H 2 = H 2 [pt(z) / Pt(0)] as 



H 2 = H 2 Pt( z ) _ u2 



10 M0) ~^ exp 

and the luminosity distance as 
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Equation (||) shows that the luminosity distance de- 
pends on a double integral over the total equation of 
state, wt{z). One integral is required to obtain the to- 
tal luminosity distance from the present back to redshift 
z. The integrand depends on H which is itself related 
to wt through the integral relation, Eq. (||). To distin- 
guish different forms of dark energy, though, we need to 
determine wq(z), the equation of state of the dark en- 
ergy component. Assuming that the Universe contains 
only pressureless matter (baryonic and cold) and dark 
energy, then wt = wqQq, where J7q is itself related to 
wq(z) through an integral relation. In particular, using 
the energy conservation analogous to Eq. (||) for the dark 
energy component alone (j?q, pq), one obtains 
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combined with Eq. (0), wt — wqVLq — wqpq/pt can be 
reexpressed as 



1 + z 

—3 J d In x (wq) 
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where £l m ,Q refers to the current values. Together with 
Eq. (0), this expression constitutes the integral relation 
between luminosity distance and wq(z) that underlies 
the degeneracy problem. 

To express the degeneracy problem quantitatively, we 
have found the maximum likelihood values of wq and w\ 
based on current SCP supernova datajj which has mea- 
sured 50 supernova out to redshift z R* 1. For simplicity, 
we have assumed wq(z) = wq + W\z; including more 
general functions of z only degrades the resolution fur- 
ther. Furthermore, we have repeated the analysis based 
on simulated data from an idealized experiment which 
measures thousands of supernovae out to redshift z — 2. 
The simulated data assumes a cosmological model with 
Q m = 0.3, J7q = 0.7 and wq = —0.7 = const. For the 
idealized experiment, the absolute magnitude measure- 
ment for a single supernova is taken to have a statistical 
variance of 0.15 and a systematic measurement error of 
0.02. The supernovae are divided into 50 bins between 
z = and z = 2 such that the statistical and systematic 
errors averaged over a bin are comparable, resulting in an 
error in cIl(z) for a given bin equal to 1.4%. We assume 
that other types of observations have constrained £l m to 
lie between 0.2 and 0.4, say, and marginalize £l m over 
that range. For the purposes of illustration, we assume 
that wq(z) is a linear function of z parameterized by wq 
and w\. As shown in Fig. lc, including more general 
forms for (z) only worsens the degeneracy problem. 

We have considered two ways of treating the system- 
atic errors. For Case I, we assume that systematic errors 
are random and uncorrelated from bin to bin and per- 
form a likelihood analysis over the 50 bins with 1.4% 
error each to determine the uncertainty in di,{z). (To 
obtain 1.4% per bin for 50 bins requires measuring thou- 
sands of supernovae.) In Case II, we assume there is 
negligible statistical uncertainty but correlated system- 
atic error of 1%. Examples of Case II errors are those 
due to calibration, dust, or evolution of supernovae. In 
this case, all models which predict d^iz) within 1% of 
the assumed cosmological model for all z between and 
2 are deemed indistinguishable. As it turns out, the 95% 
confidence contours in Case I are roughly equivalent to 
the indistinguishability region of Case II, so both cases 
give comparable results. 

For the current data, likelihood analyses based on the 
assumption that wq is constant have reported a resolu- 
tion of -1 < wq < -0.6 at the 95% leveLO When we 
repeat the analysis assuming a linear form for wq(z) and 
making no prior assumptions about wq and u>i, we find 



that neither parameter is well-determined. The degen- 
eracy obliterates the resolution of both quantities: wq 
can vary between -3.2 and -0.4 (95% confidence) and Wi 
can vary between -11.8 and 11.0. Notice the enormous 
range of wo; the 99% confidence contour includes positive 
values, so one cannot even be sure that the Universe is 
accelerating today. In cases where the Universe is not ac- 
celerating today, we can still conclude that it must have 
been accelerating recently because w\ is highly negative 
whenever wo is positive. One could argue that allowing 
large values of w\ so that wq{z) becomes much less than 
-1 or greater than +1 between z — and z — 1 (the range 
of current observations) is unphysical based on the pos- 
itivity and stability conditions that apply to most (but 
not all) forms of dark energy. Adding this theoretical 
constraint, the 95% range for wq is found to lie between 
-0.5 and -1.0 (95% confidence), in which case one may 
conclude that the Universe is accelerating today. How- 
ever, one should beware that our estimates are optimistic 
in assuming that wq(z) has only a linear and constant 
term. A safer assessment would be that, assuming posi- 
tivity and stability but no other prior about wq , one can 
conclude from present data that the Universe is acceler- 
ating today, but wq is very poorly resolved, and dwQ jdz 
can range anywhere within the imposed positivity and 
stability constraints. 

For the idealized experiment, the likelihood contours 
span a substantial range of (u>o, wi), as shown in Fig. 2. 
In the shaded ellipses the figure shows likelihood region 
if one assumes prior knowledge that Q m = 0.3 precisely. 
The contours stretch along a curve in the wo — w\ plane 
which corresponds to a near-degeneracy. It is this de- 
generacy that dashes hopes of using luminosity distance 
to measure both the current value and time derivative of 
wq. Marginalizing over £l rn expands the contours along a 
direction nearly orthogonal to the degeneracy curve. See 
the large black contours in Fig. 2. Within the 95% confi- 
dence region, wq spans a range equal to more than 35% 
of the assumed value (wq ~ —0.7), and w\ — dw/dzg 
ranges between +0.3 and -1.1. Expanding the fit to in- 
clude non linear wq(z) would expand the region even 
more (see Fig. 1). 

Another approach for measuring the time evolution of 
wq(z) is object counts, where the objects might be galax- 
ies, clusters, or halos: 4jy = n c r 2 dr, where dN is the 
number of objects in a comoving volume element r 2 drdfl 
for coordinate distance r and solid angle f2. One assumes 
that the number density of objects per comoving volume, 
n c , is constant or some known function of z. The distance 
r is related to the luminosity distance by d^^z) — (l + z)r 
(if we normalize the Friedmann-Robertson- Walker scale 
to be unity today). Hence, we can write 
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where prime represents derivative with respect to z. The 
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novel feature here is d' L (z), which entails one less integral 
of wq than d^, and, hence, perhaps an improved resolu- 
tion. Newman and DavigI3 suggest that near-future ob- 
servations of the number of dark matter halos as a func- 
tion of their circular velocity and redshift can determine 
Jldn t° within a few percent. Assuming a resolution of 
2.5% between z — 0.7 and z = 2 (more optimistic than 
their example), we find combinations of wq and wi for 
which dN/dzdQ. is indistinguishable from the assumed 
cosmological model. The indistinguishability region co- 
incides approximately with the 95% confidence contours 
in Fig. 2. Hence, objects counts are subject to essentially 
the same degeneracy problem as supernova searches. 




FIG. 2. 68% and 95% confidence contours in the 
(m?o,Wi) = {wq(z = 0), duiQ/dzo) plane for an idealized ex- 
periment which measures thousands of supernovae between 
z — and 2 = 2. The supernovae are divided into 50 bins 
with a net error of 1.4% per bin. The example assumes a 
model with Q m = 0.3, Qq = 0.7, wq = —0.7 = const., in- 
dicated by the circle. The thin shaded ellipses are 68% and 
95% confidence contours if one assumes fl m is fixed to be pre- 
cisely 0.3. The broader black contours are the result if Q m is 
marginalized over the range 0.2 to 0.4. 

Our analysis has shown that the luminosity distance- 
redshift relation and similar classical cosmological mea- 
sures are limited in their ability to resolve wq(z), even 
under optimistic assumptions (an absolutely flat Uni- 
verse, rather stringent priors for f2 m , exceptional accu- 
racy in determining luminosity distance, etc.) These con- 
clusions hold unless the errors can be reduced by at least 
1 order of magnitude or some complementary experiment 
can break the degeneracy. If a method could be found to 
reduce considerably the uncertainty in Q m - in Fig. 2, 
we assumed ft m to be in the range 0.2 to 0.4 - the de- 
generacy region shrinks along one direction. This would 
improve the resolution, but only modestly because there 
remains the second degeneracy direction shown in the fig- 
ure. It should be noted that reducing the uncertainty in 
Q m will be difficult. Most methods for determining Q m 



are dependent on some assumption about wq . In the case 
of the CMB anisotropy, for example, a degeneracy arises 
such that, for the same high-precision data, one can de- 
rive different values of Q m .depending on what assump- 
tion is made about wq(z)x3 Of course, if one assumes 
\dwQ/dz\ <C 1 (or some other prior), then deep super- 
nova searches and galaxy counts can resolve wq{z = 0) 
with impressive precision, but the value depends strongly 
on the particular theoretical assumption. 

Our conclusions also undermine claims that the super- 
nova and object count searches can determine the future 
fate of the Universe. Since the observations cannot dis- 
tinguish whether dwq /dzQ is positive or negative, they 
cannot distinguish whether wq will remain negative or 
become positive in the future, and, hence, whether the 
Universe will accelerate ever faster or cease accelerating 
altogether. 
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